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Abstract 

Here we present a study of stochastic resonance in an extended FitzHugh-Nagumo 
system with a field dependent activator diffusion. We show that the system response 
(here measured through the output signal-to-noise ratio) is enhanced due to the 
particular form of the non-homogeneous coupling. Such a result supports previous 
ones obtained in a simpler scalar reaction-diffusion system and shows that such an 
enhancement, induced by the field dependent diffusion -or selective coupling-, is a 
robust phenomenon. 
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1 Introduction 



Stochastic resonance (SR) is one of the most interesting noise-induced phe- 
nomena that arises from the interplay between deterministic and random dy- 
namics in a nonlinear system A large number of examples showing SR 
occur in extended systems: for example, diverse experiments were carried out 
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to explore the role of SR in sensory and other biological functions or in 
chemical systems (jsl). These, together with the possible technological applica- 
tions, motivated many recent studies showing the possibility of achieving an 
enhancement of the system response by means of the coupling of several units 
in what conforms an extended medium (0; S 0) ■ 



In previous works (|5|; ly) we have studied the stochastic resonant phenomenon 
in extended systems, when transitions between two different spatial patterns 
occurs, exploiting the concept of the non- equilibrium potential (NEP) 
a Lyapunov functional of the associated deterministic system that, for non- 
equilibrium systems, plays a role similar to that of a thermodynamic potential 
in equilibrium thermodynamics. Such NEP characterizes the global properties 
of the dynamics: attractors, relative (or nonlinear) stability of these attrac- 
tors, height of the barriers separating attraction basins and, in addition, al- 
lowing us to evaluate the transition rates among the different attractors. In 
another work (0) we have also shown that, for a scalar reaction-diffusion sys- 
tem with a density- dependent diffusion and a known form of the NEP, the 
non-homogeneous spatial coupling changes the effective dynamics of the sys- 
tem and contributes to enhance the SR phenomenon. 

Here we report on a study of SR in an extended system: an array of FitzHugh- 



Nagumo (jllf ) units, with a density-dependent (diffusive-like) coupling. The 
NEP for this system was found within the excitable regime and for particular 
values of the coupling strength (0). In the general case, however, the form 
of the NEP has not been found yet. Nevertheless, the idea of the existence 
of such a NEP is always underlying our study. Hence, we have resorted to 
an study based on numerical simulations, analyzing the influence of different 
parameters on the system response. The results show that the enhancement 
of the signal-to- noise ratio found for a scalar system (0) is robust, and that 
the indicated non-homogeneous coupling could clearly contribute to enhance 
the SR phenomenon in more general situations. 



2 Theoretical Framework 



2.1 The Model 



For the sake of concreteness, we consider a simplified version of the FitzHugh- 
Nagumo (0; H; lllh model. This model has been useful for gaining qualitative 
insight into the excitable and oscillatory dynamics in neural and chemical 
systems dl^ . It consist of two variables, in one hand m, a (fast) activator 
field that in the case of neural systems represents the voltage variable, while 
in chemical systems represents a concentration of a self-catalytic species. On 
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Fig. 1. We show the stable patterns that arise in the system. There is one stable 
pattern that is identically zero, i.e. Pq{x) = Pq{x) = and another which is 
non-zero Pf{x)). The patterns for the fields u{x) and v{x) are plotted in 

dashed and solid lines, respectively. The parameters, are Du = 0.3, Dy = 1, h = 2. 

the other hand v, the inhibitor field, associated with the concentration of 
potassium ions in the medium (within a neural context), that inhibits the 
generation of the u species (in a chemical reaction). Instead of considering the 
usual cubic like nonlinear form, we use a piece-wise linear version 
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+ f{u)-v + i{x,t) (1) 
+ I3u-av, (2) 



where f{u) = —u + Q{u — 0c), and C,{x,t) is a 5-correlated white Gaussian 
noise, that is {^{x,t)) = and {^{x,t)^{x' ,t')) = 2'y6{x — x')5{t — t'). Here 7 
indicates the noise intensity and 0c is the "discontinuity" point, at which the 
piece- wise linearized function f{u) presents a jump. In what follows, the pa- 
rameters a and j3 are fixed as a = 0.3 and (3 = 0.4. Finally, e is the parameter 
that indicates the time-scale ratio between activator and inhibitor variables, 
and is set as e = 0.03. We consider Dirichlet boundary conditions at a; = ±L. 
Although the results are qualitatively the same as those that could appear 
considering the usual FitzHugh-Nagumo equations, this simplified version al- 
lows us to compare directly with the previous analytical results for this system 



As in (j^, we assume that the diffusion coefficient Du{u) is not constant, but 
depends on the field u according to D^iu) = Du [1 + hQ{u — 0c)]. This form 
implies that the value of Du{u) depends "selectively" on whether the field u 
fulfills u > (hr or u < (hn- D„ is the value of the diffusion constant without such 
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"selective" term, and h indicates the size of the difference between the diffusion 
constants in both regions (clearly, ii h = then Duiu) = constant). D^iv) 
is the diffusion for the inhibitor f , that here we assume to be homogeneously 
constant. 

It is worthwhile noting that when the parameter h is large enough, under some 
circumstances the coupling term might become negative. This is what is known 
as "inhibitory coupling" (|l2l). This is a very interesting kind of coupling that 
has attracted much attention in the last years, both in neural and chemical 
context, that we will not discuss here. 

This system is known to exhibit two stable stationary patterns. One of them 
is u{x) = 0, v{x) = 0, while the other is one with non-zero values and can be 
seen in Fig. 1. We will denote with Pq'^{x) and Pi'^{x), the patterns for u 
and V fields. Further, we consider that an external, periodic, signal enters into 
the system through the value of the threshold 0c, 

<pc{t) = (pc + ^(p cos(u;t), (3) 



where ui is the signal frequency, and 6(f) its intensity. 

All the results shown in this paper were obtained through numerical simula- 
tions of the system. The second order spatially discrete version of the system 
indicated in Eqs. (1, 2) reads 

Ui = Du,i{ui-i + Ui+i - 2ui) + {Du,i+i - Du,i-i){ui+i + Ui-i) 

+f{ui)-v, + Ut) (4) 

Vi = D^{vi^i + Vi+i - 2vi) + (3ui-avi. (5) 

We have performed extensive numerical simulations of this set of equations 
exploiting the Heun's algorithm ([lih. 



2.2 Response's Measures 



Since the discovery of the stochastic resonance phenomenon, several differ- 
ent forms of characterizing it have been introduced in the literature. Some 
examples are: (i) output signal-to-noise ratio (SNR) ll^ . (ii) the spec- 
tral amplification factor (SAF) ([ist [l6h. (iii) the residence time distribution 
([itI: Eih. and, more recently, (iv) information theory based tools 



Along this paper, we will use the output SNR at the driving frequency lo. 

In this spatially-extended system, there are different ways of measuring the 
overall system response to the external signal. In particular, we evaluated the 
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output SNR in three different ways (the units being given in dB) 



• SNR for the element N/ 4 of the chain evaluated over the dynamical evolu- 
tion of MAr/4, that we call SNRi. 

• SNR for the middle element of the chain evaluated over the dynamical 
evolution ofuN/2, that we call SNR2. Having Dirichlet boundary conditions, 
the local response depends on the distance to the boundaries. 

• In order to measure the overall response of the system to the external signal, 
we computed the SNR as follows: We digitized the system dynamics to a 
dichotomic process s{t): At time t the system has an associated value of 
s(t) = 1 (0) if the Hilbert distance to pattern 1 (0) is lower than to the 
other pattern. Stated in mathematical terms, we computed the distance 
X'2['5 ■] defined by 



in the Hilbert space of the real- valued functions in the interval [—L, L], i.e. 
£2- At time t, a digitized process is computed by means of 



We call this measure SNR 



3 Results 

As indicated above, Eqs. (4) and (5) have been integrated by means of the 
Heun method (ji^. We have fixed the parameters e = 0.03, 0c = 0.52 and 
adopted an integration step of At = 10~^. For the signal frequency we adopted 
u = 271/3.2 = 1.9634295 . . .. The simulation was repeated 250 times for each 
parameter set, and the SNR was computed by recourse of the average power 
spectral density. 

Figure 2 depicts the results for the different SNR's measures we have previ- 
ously defined as function of the noise intensity 7. We adopted the following 
values: S(j) = 0.4, = 1. and N = 51. In all three cases it is apparent that 
there is an enhancement of the response for h > 0, when compared with the 
h = case, while for h < the response is smaller. 

In Fig. 3 we show the same three response's measures, but now as a function 
of h. We have plotted the maximum of each SNR curve, for three different 





ifV2 [P^{x),u{x,t)] < V2 [P^{x),u{x,t)] 
if V2 [Prix),u{x, t)] > V2 [P^ix),uix, t)] 



(6) 
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Fig. 2. SNR vs. 7, the noise intensity, for the three different measures we use. 
The parameters are 64> = 0.4, Dy = 1., uj = 27r/3.2, A'^ = 51. The different curves 
represent different values of h, showing an enhancement of the response to the 
external signal for 7 > 0. In particular it is shown: h = —2 (+), h = —1 (A), h = 
(o), /i = 1 (□) and /i = 2 (O). 

values of the noise intensity, and for 50 = 0.4, = 1., 7 = 0.01,0.1,0.3, 
Du = 0.3, and N = 51. It is clear that there exists an optimal value of 7 such 
that, for such a value, the phenomenon is stronger (that is, the response is 
larger). It is apparent the rapid fall in the response for h < 0. 

In figure 4 we show the dependance of SNR on h, for different values of the 
diffusion which depends on the activator density Du- It is apparent that the 
response becomes larger when the value of is larger. However, as was 
discussed in (j^ [gl) , it is clear that for still larger values of Du, the symmetry of 
the underlying potential (that is the relative stability between the attractors) 
is broken and the response finally falls-down. 

Figure 5 shows the results of the SNR, but now as function of D^, the activator 
diffusivity, for different values of 7, and for 50 = 0.4, = 1. and = 51. It 
can be seen that, independently from the coupling strength D^, the response 
to the external signal grows with the selectiveness of the coupling, showing 
the robustness of the phenomenon. 

Next, in figure 6, we present the results for the SNR as function of D^, the 
activator diffusivity, for different values of 7, and for 50 = 0.4, = 0.3, and 
N = 51. We see that ioT h > the response is more or less flat, however, it 
is again apparent the SNR's enhancement for /i > 0. For h < we see that 
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Fig. 3. SNR vs. h, the selectiveness of coupling, for the three dif- 
ferent measures we use. The parameters are 5(j) = 0.4, = 1., 
7 = 0.032 (O), 0.32 (□), 0.6 (o), 1.2 (□), 3.2 (x), D„ = 0.3, = 51. 
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Fig. 4. SNR vs. h, the selectiveness of coupling, for different values of Du- The 
parameters are 6(j) = 0.4, Dy = 1., 7 = 0.32, Du = 0.0 (O)>0.1 (□),0.3 (o), N = 51. 
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Fig. 5. SNR vs. -D^j, the diffusiveness in activator variable n, for the three 



different measures we use. The parameters are 



0.4, 



7 = 0.1 (O))0-32 (□),1.0 (o), while the white symbols represent h 
black ones, h = 0. The system size is = 51. 



= 1-, 
2 and the 



the system's response decays very fast with increasing D^. This effect could 
be associated to the fact (as found in those cases where the NEP is known 
(j3 0)) that in the underlying NEP the bistability is lost as a consequence of 
the disappearance of some of the attractors 



Finally, in figure 7 we depict the same three SNR's measures but as a function 
of A^, the system size. For the two measures SNRi and SNR2, we see that, for 
different values of h and 7, the response is very flat, and do not seems to be 
too much dependent on A^. It is clear that there is an increase of the response 
when h increases. At variance, for SNRp, the dependence to A^ is apparent: 
the SNR decays to zero, in a fast or slow way, depending of /i = or /i > 0. 



Here 



0.4, = 0.3, = 1., 7 



0.01,0.1,0.3. 



4 Conclusions 



We have analyzed a simplified version of the FitzHugh-Nagumo model 



111 ), where the activator's diffusion is density-dependent. Such a system, when 



both diffusions are constant (that is: Du > and = 0), has a known form 
of the NEP (0). However, in the general case we have not been able to find 
the form of the NEP (but the idea of such a NEP is always underlying our 
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Fig. 6. SNR vs. Dy, the diffusiveness in the inhibitor variable v. The parameters are 
S4> = 0.4, Du = 0.3, uj = 2 7r/3.2, = 51. 7 = 0.1 (O),0.32 (□),1.0 (o), while the 
white symbols represent h = 2 and the black ones, /i = 0. 
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Fig. 7. SNR vs. AT, the system size, for the three different measures we 
use. The parameters are Scf) = 0.4, = 0.3, = 1., 7 = 0.32, 

h = -2 {O),0 (grey squares), 2 (♦). 
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analysis) and we have to resort to an analysis based on numerical simulations. 

Through the numerical approach we have studied the influence of the differ- 
ent parameters on the system response. From the results it is apparent the en- 
hancement of the output SNR as h, the selectivity parameter, is increased. This 
is seen through three different ways of characterizing the system's response. 
We can conclude that the phenomenon of enhancement of the SNR, due to 
a selectivity in the coupling, initially found for a scalar system (j^ is robust, 
and that the indicated nonhomogeneous coupling could clearly contribute to 
enhance the SR phenomenon in very general systems. This phenomenon is 
also robust to variations of the parameter that controls the selectiveness of 
the coupling, up to a point that even in the case of inhibitory coupling the 
phenomenon holds. 

An aspect worth to be studied in detail is the dependence of the SNR on A^, 
the number of coupled units. In this way we could analyze the dependence of 
the so called system size stochastic resonance 

The thorough study of this problem will be the subject of further work. 



24) on the selective coupling. 
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